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^ , In this foUowup to Phys. Rev. D 75, 053001 (2007) [arXiv:hep-ph /0608060 we report updated 

.^1^ . constraints on neutrino mass-mixing parameters, in light of recent neutrino oscillation data (Kam- 

^ ■ LAND, SNO, and MINOS) and cosmological observations (WMAP 5-year and other data). We 

00 I discuss their interplay with the final 0u2(3 decay results in ''^Ge claimed by part of the Heidelberg- 

Moscow Collaboration, using recent evaluations of the corresponding nuclear matrix elements, and 
their uncertainties. We also comment on the 0v2l3 hmits in ^'^"Te recently set by Cuoricino, and on 
^ prospective limits or signals from the KATRIN experiment. 

PACS numbers: 14.60.Pq, 23.40.-s, 95.35.-|-d, 98.80.-k 

, Introduction. This paper is meant as a foUowup to the article [l| where, building upon previous work f3|, we presented 
constraints on the neutrino mass-squared differences (5rn^, Am^) and mixing angles (sin^0i2, sin^023, sin^^is), as 
well as on three observables sensitive to absolute v masses: the effective mass m_a in single beta decay, the effective 
Majorana mass m/3/3 in neutrinoless double beta (Oj/2/3) decay, and the sum of v masses E in cosmology — see [H, 0, 13| 
f — ' for notation and conventions. We update the results of [l| by including several new experimental inputs, largely 
T— I , presented or discussed at the recent Neutrino 2008 Conference [1]. 

in 

Neutrino oscillation updates. The Kamioka Liquid Scintillator Anti-Neutrino Detector (KamLAND) Collaboration 
has presented reactor Ve disappearance and geo-i^ results for an exposure of 2.881 kTy f5|, a factor ~4 higher than 
CD , the one we used in [ij . Following ^ , the KamLAND spectrum analysis in [l], has been upgraded [S] to include the 
00 . rates of geo-v events from U and Th decay as low-energy nuisance parameters. 

' Results from the third phase of the Sudbury Neutrino Observatory (SNO-III) recently presented at Neutrino 
^ , 2008 have been included [6] in the form of two new integral determinations of the charged- and neutral-current 
' event rates Q. Other solar v updates, with a minor impact in the global parameter estimate, include the latest 
rN \ Borexino results Q and reevaluated GALLEX data — see also p^ . 

■ The Main Injector Neutrino Oscillation Search (MINOS) Collaboration has presented accelerator disappearance 
data from 3.36 x 10^° protons on target jTl|, a factor of ~2.6 larger than previously used in In the ofhcial MINOS 
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data analysis [ll| . for any given energy profile of the i/^ survival probability P^^(-Ei/), a "beam matrix" method is 
used to map the energy spectrum from near to far, and an independent near-far extrapolation method is used as 
a cross-check [T5|. This approach can be fully implemented only within the Collaboration. For our purposes, we 
analyze the 18-bin energy spectrum ratio [Tl| by folding the function P^^{E^, Ato^, sin^ 623, sin^ ^13), with empirical 
energy resolution profiles, which mimic the near-far energy spectrum mapping of [12]. Normalization and energy scale 
systematics are treated as nuisance parameters. 

In the limit ^13 — > 0, our effective 21/ parameter fits reproduce very well the ofhcial ones as obtained by the 
KamLAND |5[, SNO-III [3], and MINOS [11| Collaborations. In our global analysis, however, we treat ^13 as a free 
parameter. 

Figure [1] displays our updated results on the mass- mixing parameters, in terms of standard deviations from the 
best fit (n„ — \/ IS.y^ after marginalization) . Table U summarizes such results in numerical form. As compared 
with [l[, the Am^ uncertainty is almost halved (by new MINOS data), and both the 8rr? and the sin^ 26*12 allowed 
ranges are reduced (by new KamLAND and SNO data). The range of sin^ 6*23 is almost unchanged. As discussed in 
[Tot , an intriguing new result is the preference for ^13 > at the level of ~1.6ct (or, equivalently, ~90% C.L.). Such 
an indication emerges from the combination of two independent hints in favor of ^13 > 0, each at the level of '^\a: an 
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TABLE I: Global oscillation analysis (2008): best- fit values and allowed ric ranges for the mass- mixing parameters. 



Parameter 


Sm^/IQ-^ eV^ 


sin^ 6i2 


sin^ 6*13 


sin^ 623 


AmVlO"^ eV^ 


Best fit 
1(7 range 
2(7 range 
3(7 range 


7.67 
7.48 - 7.83 
7.31 - 8.01 
7.14 - 8.19 


0.312 
0.294 - 0.331 
0.278 - 0.352 
0.263 - 0.375 


0.016 
0.006 - 0.026 

< 0.036 

< 0.046 


0.466 
0.408 - 0.539 
0.366 - 0.602 
0.331 - 0.644 


2.39 
2.31 - 2.50 
2.19 - 2.66 
2.06 - 2.81 



older one, found in the atmospheric v data analysis of [3[, and a newer one, coming from the small difference between 
the best-fit values of sin 26*12 in KamLAND H and SNO — a difference which is reduced for sin few %. 

Hereafter, as in [J, we shall show results at a conservative 2a (95%) confidence level, in which case only an upper 
bound can be placed on 613. 

Figure [2] shows the 2a bounds implied by the above v oscillation parameter constraints (for normal or inverted 
hierarchy) in the three planes charted by any two among the three observables (m^, mpfj, S). A measurement of any 
such quantity, coupled with the bounds in Fig. [21 provides "predictions" for the other two quantities [H, [3, S] ■ 

Cosmology updates. Within the standard cosmological model, the 5-year data recently released by Wilkinson 
Microwave Anisotropy Probe (WMAP 5y) [13, [3] constrain, by themselves, the sum of the v masses E below 1.3 eV 
at 95% C.L. This limit can be strengthened in the sub-eV range by adding further cosmological data; for instance, 
the WMAP collaboration finds S < 0.61 eV by adding Baryonic Acoustic Oscillation (BAO) and type-la supernova 
(SN-Ia) data [H. 

We consider five representative combinations of cosmological data, which lead to increasingly stronger upper limits 
on E: (1) Cosmic Microwave Background (CMB) anisotropy data from: WMAP 5y fl4], Arcminute Cosmology 
Bolometer Array Receiver (ACBAR) fis'l. Very Small Array (VSA) [l^, Cosmic Background Imager (CBI) 17] and 
BOOMERANG [l^l experiments; (2) the above CMB results plus the large-scale structure (LSS) information on 
galaxy clustering coming from the Luminous Red Galaxies Sloan Digital Sky Survey (SDSS) [I^l; (3) the above CMB 
results plus the Hubble Space Telescope (HST) prior on the value of the reduced Hubble constant, h = 0.72 ± 0.07 
[20| . and the luminosity distance SN-Ia data of [2l|; (4) the data in (3) plus the BAO data from (5) the data in 
(4), plus the small scale primordial spectrum from Lyman-alpha (Lya) forest clouds [23ll24|. 

Adopting the same procedure described in [H, based on the pubhcly available COSMOMC code [2^, we find the 
upper limits on S summarized in Table HIl and in Fig. [31 The bound in Table HIl for case (1) is dominated by WMAP-5y 
data. The results for cases (1) and (4) are in agreement with similar constraints presented in 13] and p^ . respectively, 
even if the datasets considered here for BAO and SN-Ia are different. In Fig. [31 the slight preference for S ^ at best 
fit for case (4) [and case (1)], also found in [ij], is not statistically significant. As in Ref. [l|, we find that Lya data 
have a strong impact on E, but their inclusion in global fits is debated due to systematics still under scrutiny. 
The upper limits from cases (l)-(4) (namely, S < 0.6 — 1.2 eV) should be considered as more conservative. Including 
LSS data would not significantly modify case 5, which is dominated by liya data. In the following, we shall focus on 
the two extreme cases (1) and (5). 

Oi^2/3 decay updates. The final analysis of part of the Heidelberg-Moscow (HM) Collaboration reports a 0i'2f3 signal 
in ^^Ge with half-life T^^^ = 2.23lg;^^ x 10^^ y {la errors) at a claimed C.L. > 6a [26]. The previously estimated 
T^/'2 as used in [H, was a factor of ~2 smaller. The claim is controversial, but the experimental sensitivity to 
the signal (if real) is no longer questioned [1^ . 

From a theoretical viewpoint, the Oz/2/3 nuclear matrix elements (NME) Cmm and uncertainties estimated via 
Quasiparticle Random Phase Approximations (QRPA) in ^] (as used in pj) have been recently revised [s^, HH, 



TABLE II: Representative cosmological data sets and corresponding 2a (95% C.L.) constraints on the sum of v masses E. 



Case Cosmological data set E (at 2a) 

1 CMB < 1.19 eV 

2 CMB + LSS < 0.71 eV 

3 CMB + HST + SN-Ia < 0.75 eV 

4 CMB + HST + SN-Ia + BAO < 0.60 eV 

5 CMB + HST + SN-Ia + BAO + Lya < 0.19 eV 
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especially to improve the so-called short-range correlations. We adopt for Cmm the central values and errors of [3l|, 
which agree with independent QRPA [s^, Hsj and shell model [sJl evaluations within ~2tT (see Fig. 11 of [3l[ and 
related comments therein). 

The effect of both the Tjy'2 and the NME updates for ^^Ge is to lower the central value — and to enlarge the errors — 
of the effective mass parameter m^^ = m1/CmmT^'/2- taking logs (in base 10) to linearize error propagation, we 
have log(T0/'^2/y) = 23.35 ± 0.16 {2a) from [H and log(C™^/y"i) = -12.82 ± 0.48 (2cr) from so that 

logimpp/eV) = -0.54 ±0.26 (HM claim, 2a) , (1) 

where the experimental error and the (dominant) theoretical error have been added in quadrature. 

The Cuoricino experiment, which does not find Oi^2/3 decay signals in ^^°Te, quotes > 3.1 x 10^** y at 90% 

C.L. [28|, or > 2-5 x 10^^ y at 95% C.L. [35|. Using the latter limit as logiT^J'^/y) > 24 A, and the ^^Oxe NME 

estimate log(C„„i/y"i) = -12.27 ± 0.28 {2a) from [3l|, we get 

log(m/3/3/eV) < [-0.63, -0.07] (Cuoricino, 2a) , (2) 

where the range due to the 2a uncertainty of the NME is explicitly reported. 
A comparison of the corresponding m/3/3 ranges {2a), 

0.16 < mpp/eY < 0.52 (HM claim) , (3) 
< mpp/cV < 0.23 (Cuoricino, "favorable" NME) , (4) 
< mpfs/eV < 0.85 (Cuoricino, "unfavorable" NME) , (5) 

shows that current Cuoricino data may or may not disfavor a fraction of the HM range for mpf^ at 2a, depending 
on the (still quite uncertain) value of the ^^°Te 0i^2/3 NME. A similar conclusion (albeit with somewhat different 
preferred ranges for TO/3/3) has been reached in [28j . Therefore, the 0^2(3 claim [26| remains an open issue at present, 
and we shall consider the possibility that it corresponds to a real signal. 

Discussion. Figure 3] shows the regions allowed at 2a in normal and inverted hierarchy (slanted bands) by the 
combination of oscillation results with the first dataset in Table HIl (CMB). in the plane spanned by (E, mpp). This 
is the most conservative case, with the weakest limits on S, and the largest overlap between the regions separately 
allowed by oscillation-|-CMB data and by the 01/2(3 claim. The results of a global fit are shown as a thick black 
wedge in the upper right part of the figure. [The combination includes the current limit m/3 < 1.8 eV {2a) [H which, 
however, provides only a minor contribution.] Such global combination would correspond to nearly degenerate masses 
in the range 

mi ~ TO2 ^ m3 e [0.15, 0.46] eV {2a) . 

In this case (degenerate spectrum), the preferred range for effective neutrino mass in f3 decay would also be 
TTip £ [0.15, 0.46] eV. In the upper half of this range, the KArlsruhe TRItium Neutrino (KATRIN) /3~ experiment 
could make a 5a discovery, according to the estimated sensitivity [365 • ^ 3ct evidence could still be found in KATRIN 
for rufj ~ 0.3 eV. Below this value, the sensitivity would be rapidly degraded, and only upper bounds could be placed 
for rriff < 0.2 eV [36ll. The possibility of reaching a ~0. 1-0.2 eV sensitivity with a different approach to /3 decay is 
being discussed (371 . 

If the cosmological dataset (1) were replaced by the datasets (2)-(4) in Table [Til the overlap region between the Oi^2/3 
band and the oscillation-|-cosmological bands in Fig. [4] would shrink (not shown), but would not disappear. Therefore, 
within the standard framework and the present uncertainties, the 0^2(3 claim clashes with oscillation-|-cosmological 
data only if the latter include Lya data. 

Figure [S] is analogous to Fig. 21 but refers to the fifth dataset in Table [Til fall cosmological data, including Lya). In 
this case, the allowed regions do not overlap and cannot be combined, since the relatively strong cosmological limit 
S < 0.19 eV implies rn/3/3 < 0.08 eV, in contradiction with Eq. ([3]). Solutions to this discrepancy would require that 
either some data or their interpretation are wrong. 

In conclusion, important pieces of information are being slowly added to the puzzle of absolute v masses. In 
this foUowup to [l|, we have discussed the most recent oscillation and nonoscillation updates in the field, after the 
recent Neutrino 2008 Conference Q. Oscillation parameters are robustly constrained, and an intriguing indication for 
6*13 > emerges, as summarized in Fig. [T] and in Table [D Concerning nonoscillation observables, despite some recent 
experimental and theoretical progress, a coherent picture remains elusive. In particular, the Oi^2/3 claim is still under 
independent experimental scrutiny, and it may be compatible (Fig. 3]) or incompatible (Fig. [5]) with the cosmological 
bounds (Table |TT| , depending on data selection (especially Lya) . A confident assessment of the v mass scale will 
require converging evidence from at least two of the three observables {rrif}, mpp, E) within the bands allowed by 
oscillation data in Fig. [2l 
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Synopsis of global 3z/ oscillation analysis 
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FIG. 1: Global oscillation analysis (2008): Bounds on the mass-mixing oscillation parameters, in terms of standard deviations 
from the best fit. Note the 1.6a preference for 613 > 0. 
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FIG. 2; Bands allowed at 2a by neutrino oscillation data, in each of the three coordinate planes of the parameter space 
(m/3, m^f), E), for both normal and inverted hierarchy. 
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FIG. 3: Cosmological constraints on the sum of neutrino masses (S). Standard deviation curves for the five datasets in Table Ull 
1 (dotted), 2 (dashed), 3 (dot-dashed), 4 (long dashed), and 5 (solid). 
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V oscilL + ^ + 0i/2^ claim + CMB 
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FIG. 4: Global combination of oscillation plus CMB data (case 1 in Table with the ^v2(i decay claim, in the plane (E, mpp). 
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FIG. 5: Bounds from oscillation plus all cosmological data (case 5 in Table . contrasted with the 0/^2/3 decay claim, in the 
plane (E, m/sff). 



